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We report the electronic properties of two-dimensional systems made of graphene nanoribbons 
which are patterned with ad-atoms in two separated regions. Due to the extra electronic confinement 
induced by the presence of the impurities, we find resonant levels, quasi-bound and impurity-induced 
localized states, which determine the transport properties of the system. Regardless of the ad-atom 
distribution in the system, we apply band-folding procedures to simple models and predict the 
energies and the spatial distribution of those impurity-induced states. We take into account two 
different scenarios: gapped graphene and the presence of randomly distributed ad-atoms in a low 
dilution regime. In both cases the defect-induced resonances are still detected. Our findings would 
encourage experimentalist to synthesize these systems and characterize their quasi-localized states 
employing, for instance, scanning tunneling spectroscopy (STS). Additionally, the resonant transport 
features could be used in electronic applications and molecular sensor devices. 


I. INTRODUCTION 


Two-dimensional nanostructures are promising candi¬ 
dates for novel application at the nanoscale. Due to its 
dimensionality, it is possible to manipulate the entire 
system through several experimental techniques. These 
nanosystems are suggested for sensing applications be¬ 
cause they have large surface-volume ratio, high elec¬ 
tronic mobility, and externally tuneable conductivity. 
One of its applications is related to their use in high- 
precision molecular and magnetic sensors, which requires 
robust states at well-defined energies.®® In 2D electron¬ 
ics, the control of the electronic and magnetic properties 
can be achieved by specific point defects, such as va¬ 
cancies, Stone-Wales defects and doping with ad-atoms, 
among others. 7 11 Previous theoretical works propose 
electron confinement by doping 2D graphene with hy¬ 
drogen atoms, which is used to manufacture narrowed 
graphene ribbons, quantum dots or junctions without the 
need for cutting or etching the system. 12 ®^ Furthermore, 
it was experimentally shown that the adsorption of hy¬ 
drogen atoms on graphene can be controlled by scanning 
tunneling microscope to decorate with patterns the sys¬ 
tem, which remains stable even at room temperature?® 
Graphene nanoribbons ideally could be today synthe¬ 
sized by unzipping carbon nanotubes. This method al¬ 
lows provides narrow graphene nanoribbons suspended 
between Au (111) contacts or deposited on a substrate 
such as SiC or BN? 6 19 Such nanoribbons would be later 
protected from hydrogen deposition by a controlled coat¬ 
ing, for instance, of PmPV (m-phenylenevinylene-co-2,5- 
dioctoxy-p-phenylenevinylene) at two different regions. 19 
The partial hydrogenation could be achieved by either di- 
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re ct exp osure or exposing the sample to a beam of atomic 
H 1M11U Then, scanning tunneling microscopy (STM) and 
scanning tunneling spectroscopy (STS) could be used to 
characterize the quasi-localized states induced by the two 
H regions. 2 ^Note that STM experiments are able to con¬ 
trol the H deposition in graphene sublatticed 23 ! In this 
way, systems as those shown in Fig. [I] can be produced. 



FIG. 1: (Color online) (a) Proposed system based on GNR 
with ad-atoms in two separated regions, (b) Schematic view 
of ad-atoms in ordered configurations. The length L (in units 
of 3 acc ) is the separation between the two doped regions. 
The I and D notations stand for individual and double hydro¬ 
gen column. Red spheres represent hydrogen ad-atoms. The 
effective system length L e // is given by the shaded region. 
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In graphene, the hydrogen-like ad-atoms are adsorbed 
on top of the carbon atoms through covalent bonds. 
The local hybridization of graphene changes from sp 2 
to sp 3 -like, so th at the low-energy conducting 7r-bands 
are affected J 24 * 25 * When hydrogen-like atoms are binding 
on top of carbon atoms, the graphene sublattice symme¬ 
try is broken, and thus, localized states appear near the 
Dirac point. These states at zero energy are well-known 
in carbon nanostructureeP 27 and are comparable, in the 
low-energy range to a vacancy in 7 r models ? 28 1 

Charge carriers in graphene have a linear energy 
dispersion^ then electron scattering obey the Klein’s 
paradoxPSI Klein showed that when the potential is of the 
order of the electron kinetic energy, the barrier for mass¬ 
less fermions becomes transparent. In practical terms, 
the Klein’s paradox implies that car riers in graphene can 
not be confined by potential barriers ? 30 31 In this context, 
the question that motivates this work is the following: 
Is it possible to induce and control resonant states, ex¬ 
tended across the whole system, by using impurities in 
graphene?. 

In this paper we study the electronic and transport 
properties of graphene ribbons patterned with hydrogen- 
like ad-atoms at two separated regions, as it is shown in 
Fig(lJ This approach can be used as a model for any 
singly covalent adsorbed molecules in graphene. We are 
describing the H interaction on top of C atoms using 
a tight-binding Hamiltonian, which are compared with 
density functional theory (DFT) calculations. We cal¬ 
culate the local density-of-states (LDOS) and the two- 
terminal conductance of the system, and we are con¬ 
sidering the effects of the H-adsorbate impurity density, 
the spatial distribution and the sublattice balance at the 
doped regions. We find not only the expected impurity- 
induced localized states, but other kinds of electronic 
states such as resonant and interface states. These levels, 
which determine the transport properties of the system, 
appear at well-defined energies that we can unravel, and 
additionally infer their spatial distribution, by applying 
band-folding models. 

To bring our results into contact with experiments, we 
examine a random distribution of hydrogen adsorbate 
atoms at low concentrations in the two doped regions. 
The average LDOS still exhibits a series of defined peaks 
at some energy values, which behave similarly to the or¬ 
dered ad-atom arrangements. We are also study a ribbon 
over a substrate, modeled by a constant staggered poten¬ 
tial (a mass term in the Dirac equation) which open a 
gap in the LDOS and conductance curves. In the latter 
system, the behavior of the electronic states is still be 
determined by applying the same kind of simple folding 
models. 

Note that our proposal goes further than the ex¬ 
isting reports, which have considered periodic bound¬ 
ary conditions and/or superlattice heterostructures, 
where the electronic properties and the transport are 
calculat ed paralle l to the confinement generated by 
barriers j 12 l 1 ^ l 32 l 3 H 


II. THEORETICAL DESCRIPTION OF THE 
SYSTEM 

The low-energy electronic properties of graphene are 
well described within a single 7r-band approximation in a 
nearest-neighbor tight-binding scheme. The Hamiltonian 
of the system can be written as 

U = -70 Y c l c m+SHdid n -'jH (4 n d n + c p ji ), (1) 

( l,m) 

where the hopping energy between nearest-neighbor car¬ 
bon atoms is 70 = 2.75 eV, Cm((4n) * s th e C annihilation 
(creation) operator on the ra-th site of graphene lattice, 
7 h is the coupling energy between the carbon and the 
hydrogen ad-atom, Eh is the hydrogen on-site energy, 
d n (dJJ is the H annihilation (creation) operator on the 
adsorbate site, and c Pn is the host carbon atom bonded 
to the hydrogen atom. We use the set of tight-binding 
parameters 7 h = 5.72 eV = 2.08 70 and Eh — 0 deter¬ 
mined from first-principles band structure calculations of 
hydrogenated graphene?®^. It is worth to mention that 
the effects of considering eh 7 ^ 0 are not being considered 
because preliminary tests show that this value split en¬ 
ergy states outside of the energy window studied in this 
paper. Details about the calculations of the conductance 
and LDOS are presented in the Appendix[A] We focus on 
armchair nanoribbons, which at the edges mix sublattices 
and valleys, so that the spin the degree of freedom entails 
to have doubly degenerate bands, in contradistinction to 
zigzag ribbons . 37 


III. RESULTS 

In this section we present results of LDOS and conduc¬ 
tance of armchair graphene nanoribbons (AGNR) in the 
presence of ordered and random distributions of hydro¬ 
gen ad-atoms in two separated regions of the ribbon. Be¬ 
sides, we explain the LDOS behavior by applying folding 
of certain bands of simple effective models. The LDOS 
is projected in all atoms of the system, marked by the 
outer dashed lines in Fig. [l] 


A. General Trends on the Local Density of States 
and Conductance 

We chose a metallic AGNR in the range of the ultra¬ 
narrow ribbons, recently synthesized ! 38 The conductance 
curves of these narrow ribbons exhibit broad plateaus at 
low energy, making easier the visualization of the dop¬ 
ing effects. Nevertheless, we would like to note that our 
findings in the LDOS and conductance curves are gen¬ 
eral and appear in wider metallic and semiconductor rib¬ 
bons. The considered separation between the two doped 
regions is L = 15 = 45 acc and the AGNR width is 








3 


N = 11 = 5\/3acc- In this nanoribbon, we have consid¬ 
ered the three ordered arrangements of adsorbed hydro¬ 
gen atoms following the scheme of Fig. [l] 

We start the analysis by examining the LDOS as a 
function of the Fermi energy, as it is shown in Fig. 
[ 2 ] Our first finding in the LDOS curves is a strong 
peak at zero energy, even for semiconductor ribbons, 
the so-called mid-gap state. This impurity-induced 
state is attributed to the repulsive short-range poten¬ 
tial around monovalent impurities (atoms or molecules), 
which ar e stro ngly bonded to one of the carbon atoms in 
graphene lilM 27 IMSil Generally, this mid-gap state could 
be related to the magnetic properties of the material, 
which it has been widely studied and is beyond the scope 
of this work.® 



FIG. 2: (Color online) Local density of states (black line) 
and conductance (red line) as a function of energy for different 
hydrogen distributions, separated by a length L — 15 on a 
metallic AGNR N — 11. 

Besides of the impurity-induced peak at E = 0, we 
observe common features in the LDOS curves for all con¬ 
figurations in Fig. [ 2 ] At energies below 0.4 70 , defined 
by the first van Hove singularity of the pristine N = 11 
AGNR, the LDOS curves present a series of well-defined 
sharp peaks homogeneously spaced. The energy separa¬ 
tion and the total number of those peaks are determined 
by an effective system length L e ff{> L), defined as the 
mean distance between the doped regions, as it is shown 
in shaded on Fig. [l] In this sense, the separation between 
I-configuration resonances is shorten than those of the D- 
configurations, since L\^ > L^ D2 \ These states get 
sharper and closer in energy as L e ff is increased, and it 
is possible to follow continuously their evolution as the 
system length is enlarged. If the ribbon width N is in¬ 
creased, we can still find resonant levels in the low energy 
range similar to the presented in Fig. [2j independently 
if the ribbons are metallic or semiconductor. For wider 


metallic ribbons, the low-energy range is minor; never¬ 
theless, is still possible to observe the same patterns of 
LDOS. For semiconductor ribbons we have observed the 
impurity mid-gap state at E = 0, and resonant levels 
above the in herent g ap of the ribbon. 

This quantum well-like behavior is, in fact, a hall¬ 
mark of the resonances ar ising from confinement between 
the hydrogenated regions S 21 MJ The re gular periodicity of 
these states allows us to follow resonances above energies 
of 0.4 70 , which are mixed with other kinds of energy lev¬ 
els. In the following sections, we will explain this LDOS 
behavior by using simple band-folding models. 

In Fig. [2] it is also presented the two-terminal conduc¬ 
tance of the system as a function of energy (red online). 
The hydrogen ad-atoms impose an extra spatial elec¬ 
tronic confinement, which produces some resonant states 
between the doped regions, that act as open channels for 
electron transmission. The conductance profiles present 
a series of well defined peaks becoming from strongly 
quasi-localized states, which acquire a finite lifetime due 
to their indirect coupling through the continuum spectra 
of the leads, in analogy of the mechanism of formation of 
Bound States in the Continuum (BICs) in high symmet¬ 
rical systems.® 

For the I-configuration, the doped regions present se¬ 
lective doping where hydrogen ad-atoms occupy one of 
the grapehene sublattices. Such arrangement may arise, 
for instance, when long-chain hydrocarbon molec ules are 
attached periodically on certain graphene sites 
This selective doping favors the hybridization between 
the resonant states and the continuum spectra of the 
leads, allowing a better conduction through the system 
in the low-energy range and even at higher energies. This 
clear resonant behavior, where the conductance curve 
presents well defined conduction energy channels, could 
be used for novel electronic applications, such as spin 
filter or molecular sensor devices. 

In the figure [2] we observe that the conductance curves 
for the D1 and D2 configurations are clearly different in 
comparison with those of I-configuration. For low energy, 
although the sharp states are also present in the LDOS 
curves, the D1 and D2 configurations do not conduce be¬ 
cause the C-H bonds in each doped region are balanced in 
both graphene sublattices (see Fig. [l]). The D1 configu¬ 
ration is a poor conductor (similar to an insulator) for all 
energies. This ad-atom distribution has pairs of contigu¬ 
ous carbon atoms with adsorbed hydrogen ad-atoms in 
different sublattices and, as a consequence, the electronic 
transport is mostly blocked. In fact, at low energies the 
effect on the conductance of the C-H bond is similar to 
an ideal vacancy since the 7r-type bonds are no longer 
available for conduction.® The D2-configuration favors 
the electron transmission at energies over 0.4 70 , which 
is defined by the bottom of the first parabolic conduction 
sub-band, as in the case of a line-defects in graphene.® In 
particular, for this configuration, the attached hydrogen 
atoms are separated by carbon atoms, therefore, localized 
states appear between the impurity positions which en- 

































































































B. Characterization of Resonances and 
Inter face-states by Band-folding Models 
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FIG. 3: (Color online) (a) Model used for folding the elec¬ 
tronic structure in the k x -direction. The atoms in the unit 
cell are highlighted, (b) Band structure of the pristine metal¬ 
lic ribbon N = 11. (c) Zoom of the low-energy bands with the 
k x -space quantization for a length L — 15, (vertical dashed 
lines). The horizontal lines correspond to the energies where 
the first band intersect the quantization lines, (d) LDOS of 
the D2 configuration for L — 15. 


hance the conductance through their hybridization with 
the resonant states of the system. 

It is noteworthy that, for all configurations, the con¬ 
ductance curves above the first van Hove singularity are 
dominated by parabolic bands, well-known in semicon¬ 
ductors, so they can be outlined using a toy model for 
non-relativistic electrons incident on a double-barrier po¬ 
tential within the Schrodinger theory.^ 


In this section, we unravel the tendencies of the LDOS 
curves shown in Fig. [2] We use the electronic structure 
of simplified models, which are folded in the reciprocal 
space according to the width and length of the system. 
We find a simple procedure to predict the energies of 
the resonances and their spatial distribution across and 
along the system. In what follows, we show results for the 
D2 distribution; nevertheless, other configurations can be 
equally explained with this approach. 

We start the analysis by explaining the sharp states 
exhibits in the LDOS curve of the considered system 
(Fig. 1 (d)), which arise from the metallic band of the 
ribbon. We have used a simplified model of a pristine 
N = 11 AGNR with the unit cell shown in Fig. [3] (a), 
and the corresponding band structure displayed in Fig. [3] 
(b) and (c). Hydrogen ad-atoms impose different bound¬ 
ary conditions at both ends of the system. Regarding 
the doping distribution per sublattice at the ends, the 
D2-case present a crystallographic spatial inversion sym¬ 
metry which determines the properties of the confined 
electron wave-function. The effective length L e ff of the 
system is defined by the distance between the centers of 
the doped regions, which define the band folding pro¬ 
jected along the ribbon axis. For each sublattice, the 
confined states are similar to standing waves in a pipe 
with open end (details in Appendix 0, and they have a 
discrete set of n allowed k x -values, given by: 


K\ = 


7r 

L eff 



( 2 ) 


where L e ff is the effective length (note that L e ff > L , 
presented as the shaded regions of Fig. [I]), and the n 
index is a positive integer greater than zero. 

By using the Born-von Karman boundary conditions, 
the intersection of the allowed fcJFvalues (vertical red on¬ 
line) with the envelope band (diagonal black online) de¬ 
termine the energies of the series of equidistant peaks in 
the LDOS curves of Fig. [3] (d). Here, the horizontal 
lines (green online) indicate the energies of those inter¬ 
sections, which are completely aligned with the LDOS 
peaks. The spatial distributions of the lower metallic 
states are shown in the Appendix [Bj which are in a good 
agreement with those calculated by density functional 
theory (DFT). In this sense, we are testing our calcu¬ 
lations against more sophisticated methods included in 
Appendix [Cj 

It is noteworthy to mention that we can apply the 
same folding procedure to semiconductor ribbons. We 
find that the k x -space discretization also accurately pre¬ 
dicts the energy at which the resonant sharp states ap¬ 
pear. The difference with the metallic ribbons is that now 
these states are no longer equidistant, due to the curva¬ 
ture of the parabolic envelope band, which depends on 
the ribbon width. Even though resonant states are pre¬ 
dicted with high accuracy, regardless of the character of 
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the doped regions). The interface states are identified 
by using an effective scheme, where each doped region of 
the system is individually modeled as an infinite line of 
defects forming a crystallographic discontinuity between 
two semi-infinite graphene sheets, as it is depicted in Fig. 
H (a). Considering the unit cell of the effective model 
(highlighted in this figure with blue online) we have ap¬ 
plied the Bloch’s theorem in the perpendicular direction 
of the system with a lattice constant d = ^/3 acc • We 
consider nearest-neighbors interactions in the vertical di¬ 
rection, with a /^-dependent hopping parameter. Thus, 
we have calculated the /^-projected density of states 
(pdos), given by: T>(E,k y ) = -^Tr [Im G c (E,k y )\, 
where Gc is the Green function defined in the Appendix 
|A| Note that the V (E,k y ) now is /^-dependent, as shown 
in Fig. [4] (b) for the D2-configuration. 

We observe the expected interface states at E = 0 in 
the range 2/3 < \k y d/7r\ < 1. These flat bands, corre¬ 
spond to localized edge states and are associated with 
the characteristic zigzag ends of semi-infinite graph ene 
sheet in the simplest non-interacting approach . 50 51 Al¬ 
though these states around E — 0 are not the main 
issue of this work, note that in DFT calculations they 
show spin polarization ^ 52 results which are the source 
of recent controversy due to the possible effects of spin 
contamination.^^ 

We also find an electron-hole symmetry of the bond¬ 
ing and anti-bonding interface bands produced by the 
impurities, which lie between the graphene continuum 
of states. At the T-point the splitting of these levels 
is approximately 70 , and it decreases as the horizontal 
distance between the C-H pairs within the barrier is in¬ 
creased. This energy splitting of the interface states only 
appears when the sublattices are locally balanced with 
the same number hydrogen pinned at each graphene sub¬ 
lattice, otherwise, the interface bands are not split, and 
the states remain at E = 0, extending from \k y d/7r\ = 0 
to 1 . 

For AGNRs, applying the folding procedure in the 
perpendicular direction generates a set of n allowed k 
values (vertical lines in Fig. |4](c)), which are given by: 


FIG. 4: (Color online) (a) Effective model used for the band 
folding in the k y -direction, (b) Projected LDOS (pdos) in the 
reciprocal space (ky -direction), (c) Zoom of the pdos with 
the k y -quantization lines (vertical lines) for a N = 11 AGNR. 
The horizontal lines correspond to the energies where the pro¬ 
jected pdos and the quantization lines intersect, (d) LDOS of 
the D2-configuration for a N = 11 AGNR and doped regions 
separated by a length L = 15. 


the bands for any system length, the precision decreases 
with the number of sub-bands. 

Although above the first van Hove singularity we can 
still identify states from the metallic envelope-band, the 
LDOS curves display other broader states which arise 
from the folding of parabolic sub-bands of higher ener¬ 
gies or from the interface states of the system (within 


K\=n^, (3) 

where n = 1 , 2 ,... and W = ^ (N — 1) acc corre¬ 
sponds to the ribbon width. The intersection between 
the obtained /^-values with interface bands determines 
the energy at which these states appear in the LDOS 
curves of the system, as it is shown in Fig. [4] (c) and 
(d). As the interface bands coexist with the graphene 
continuum states, for k y > 0.4, they become broader 
indicating shorter lifetimes. Due to the hybridization be¬ 
tween these states (Fig. [4] (c)), the electron conduction 
at these energy levels is allowed, bottom panel of Fig. [2j 
The spatial distribution of interface states is presented in 
the Appendix [B] These interface states appear in other 
carbon-based nanostructures such as grain boundaries in 
graphene and jucntions in carbon nanotubesP* 58 It is 
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worth mentioning that these interface states, usually con¬ 
sidered as a weakness in device performance, may provide 
a method for achieving diode behavior at the nanoscale.^ 

C. Dependence on the Ribbon Width 

We focus on the effects of the AGNR width on the res¬ 
onant and interface-states. We choose a system length 
L = 15 for which we have applied the folding procedure 
in the /^-direction in Fig. [ 5 ] (a). For metallic AGNRs, 
energy states stemming from the linear metallic band are 
shown in black marks in this panel. These states are reg¬ 
ularly distributed in energy determined by the system 
length L. For the parabolic sub-bands, the allowed states 
are given by red marks, and we identify them with some 
peaks of LDOS. In the upper-right panel, we indicate 
with blue bars the energy of the states arising from the 
metallic envelope-band and with red bars those states 
emerging from the folding of the parabolic sub-bands. 
Note that for an AGNRs width of N = 11, shadowed in 
left panels, there is a good agreement between the LDOS 
and the metallic resonances obtained by the folding pro¬ 
cedure. Besides, states coming from parabolic bands are 
more difficult to follow since they are spread in the doped 
regions and, as a consequence, their energies are shifted. 



FIG. 5: (Color online) (a) Energy position of the intersection 
of the bands and the k x -values as a function of the AGNR 
width, for a system length L — 15 (b) Interface states arising 
from the £ 7 -folding for any AGNR width. As a comparison, 
we have included the LDOS and conductance for the width 
N — 11 and a system length L — 15 (shaded region). 

By applying the folding procedure in the /^-direction, 


we identify the energies of the interface states as a func¬ 
tion of the AGNR width, as shown in Fig. [ 5 ] (b). These 
interface states are observed for energies greater than 
O .570 independent of the system width. Above 0.6 70 
in the energy range of the hybridization between inter¬ 
face states and continuum bands of graphene, the inter¬ 
face peaks become broader and contribute to the conduc¬ 
tance. In the bottom-right panel, we compare the LDOS 
and conductance with the blue marks belonging to the 
k y band-folding for an AGNR of width TV = 11. 

D. Comparison with Experimental Conditions 

In this section we aim to get into contact with the 
experiments; therefore, we focus on how these quasi- 
localized, resonant and interface states behave in gapped 
graphene or under random distributions of adsorbed hy¬ 
drogen at low concentration. 



FIG. 6 : (Color online) (a) Low energy bands for gapped 
graphene (red online), the h™ quantization for the effective 
length L = 15, has been included as vertical dashed lines. 
The horizontal lines determine the energies where the first 
band of the gapped case and the quantization lines intersect, 
(b) the density of states and (c) two terminal electronic con¬ 
ductance of the D2 configuration for a separation L — 15. For 
comparison, the gapless case has been included by black lines 
in all the panels. 


1. Gapped Graphene 

In previous sections, we have focused on the electronic 
properties of systems based on nanoribbons made of 
gapless graphene. However, most applications in elec¬ 
tronics require a band gap, so a some actual devices 
are manufactured over gapped graphene. The interac¬ 
tion of the graphene layer with the substrates, such as 
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hexagonal boron nitride (hBN) or silicon carbide (SiC), 
opens band gaps, which depend on the lattice mismatch 
and the relative angle between the graphene and sub¬ 
strate lattices E33EID Additionally, two-dimensional ma¬ 
terials with hexagonal symmetry such as hBN, silicene 
or germanium can be described by using a graphene- 
like tight-binding Hamiltonian with an additional mass 
t erm E3H66]i n thi s context, we model the gapped graphene 
by considering a constant staggered potential, namely, a 
local potential that acts with the opposite sign in the two 
sublattices. In this model, the mass of the electrons is 
a parameter corresponding to the difference | Va — Vb\\ 27 ^ 
which in the case of graphene over a SiC substrate cor¬ 
respond to Va(b) = ±0.l7o.® 

The electronic states of the system are modeled by 
using the same unit cell described in Fig. i (a). By ap- 
plying the folding procedure along the k x -direction it is 
possible to calculate a set of allowed discrete values k™ 
defined in Eq. ©• Note that these discrete momentum 
values are the same of the gapless graphene system, be¬ 
cause the mass term is diagonal, and does not affect the 
translational symmetry of the system. In Fig. [6] (a), for 
metallic N = 11 AGNR, we compare the energy bands of 
gapless and gapped graphene, which are denoted by the 
black and red lines respectively. The main effect of the 
mass term is to transform the previously linear metal¬ 
lic band into the first parabolic band, which defined the 
gap of the system. Near the Fermi level, the energy dis¬ 
tance between the LDOS peaks, determined by the fold¬ 
ing procedure, changes because the envelope band is no 
longer linear, as seen in Fig. [6] (b). At higher energies, 
the separation between the peaks of the LDOS gets close 
to previous values. The conductance curves are slightly 
affected since the effects of the mass term are mostly ob¬ 
servable near the Fermi level, Fig. [6] (c). The impurity 
induced mid-gap state previously located at E = 0 is now 
split in to series of peaks within the gap defined by the 
staggered potential. These states preserve the electron- 
hole symmetry only when the doped sublattice balance is 
conserved. We find that these defect states present sub¬ 
lattice localization in the doped region, in the vicinity 
around the hydrogen ad-atoms. All the main character¬ 
istics of these in-gap peaks, such as energies, degeneracy 
and spatial distribution for gapped graphene, can be in¬ 
ferred in detail by using band-folding analysis similar to 
the presented before in Fig. [4j 


2. Randomly Adsorbed Hydrogen at Low Concentrations 

We wonder about how robust is the effect of confine¬ 
ment in the presence of a random distribution of H im¬ 
purities. We consider low concentrations of H ad-atoms 
randomly distributed in the two doped regions. The per¬ 
centage of doping is defined as the ratio (Nh/Nc) x 100, 
where Nh is the number of H ad-atoms and Nc is the 
total number of C atoms in the doped region. Finally, 
we impose that the number of hydrogen ad-atoms is the 


same in each doped region. 

Figure [7] shows the conductance and LDOS for the N = 
11 and L = 15 AGNR as a function of the Fermi energy, 
averaged over 100 independent random realizations for 
low coverage. The concentration of random hydrogen 
doping strongly affects the electronic properties. As the 
concentration increases, the conductance of the systems 
rapidly decreases due to the disorder generated by the 
random adsorption at the doped regions. 

We note that at low energy and low concentration of 
H-atoms, there are well-defined energy levels in averaged 
LDOS and conductance curves, shown in Fig. [ 7 ] (a) and 
(b). These states resemble to those belonging to the 
ordered configurations; however, the energy separation 
between them is the half separation of the previous re¬ 
sults, which is related with the statistical correlation of 
the considered disorder. We find that the number of res¬ 
onance peaks increases with the average separation be¬ 
tween doped regions. Since the low-energy oscillations 
appeared both in conductance and LDOS, it seems that 
scanning tunneling spectroscopy could be used to detect 
those oscillations, as they were already measured in car¬ 
bon nanotubes pH 

Additionally, regardless of the distance between the 
doped regions, we observe the Fano-like resonance d 45 * 67 * 
in the average conductance, at energies corresponding to 
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FIG. 7: (Color online) (a) Average conductance and (b) av¬ 
erage LDOS as a function of the energy for a iV = 11 AGNR 
separated by L = 15 for several hydrogen coverages. Black 
lines indicate H coverage per doped region of 2.3 %; red, 4.5 
%; blue, 9.1 %. The conductance of the pristine case is shad¬ 
owed. The inset of the (a) panel includes the normalized 
conductance at the first van Hove singularity. 













the van Hove singularities of the pristine system. This 
effect could be used to determine indirectly the width of 
the ribbons employed in the system. 

IV. FINAL REMARKS 

In summary, we have shown that is possible to con¬ 
fine electrons in a GNR-based system. We accomplish 
that, by theoretically studying hydrogen-like ad-atoms in 
two separate regions. Our results show different types of 
electronics states, such as resonant and interface states, 
which determine the transport properties of the system. 
We have found evidence of the induced quantum well-like 
behavior in case of gapped graphene or regardless of the 
hydrogen distribution in the doped regions, whether they 
are highly ordered or diluted random distribution. 

We are able to predict the energies and the spatial dis¬ 
tribution of these electronic states by applying a simple 
band-folding scheme to some basic models. For the reso¬ 
nant states, we have apply the Born-von Karmarn bound¬ 
ary conditions to the metallic envelope-bands of the pris¬ 
tine AGNR. Thus, the energies of the resonant states are 
equally spaced. Furthermore, this band-folding scheme 
applies to the parabolic envelope band in both gapped 


graphene and semiconductor ribbons. Perpendicular to 
the doped ribbon, the folding scheme drives to interface 
states separated in two groups: states that hybridize or 
not with the continuum band. We found that hybridized 
states contribute with broad peaks to the electronic con¬ 
ductance. These proposed band-folding was successfully 
compared with density functional theory calculations. 

Our results would encourage experimentalist to make 
these kind of systems and characterize their quasi- 
localized states, for instance, employing scanning tun¬ 
neling spectroscopy. The obtained resonant electronic 
transport behavior could be externally manipulated by 
applying gates and/or side potentials, which allows to 
control the conductance of the system. This fact can be 
used in novel electronic applications and molecular sensor 
devices. 
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Appendix A: Calculation of Conductance 

We include the technical details concerning to the cal¬ 
culation of the LDOS and the two terminal conductance 
using the Landauer formula and the surface Green’s func¬ 
tion matching formalism.® ^ This method divides the 
system in three blocks. There are two semi-infinite leads 
made of AGNR, left L and right R described with Hf and 
Hr Hamiltonian respectively. The finite central part is 
described by the Hamiltonian matrix He , which includes 
the two doped regions and the separation between them, 
as depicted in Fig. [l] Thus, the total Hamiltonian is 
given by: 

H = H c + H r + H l + V lc + V rc , (Al) 

where Vlc ? Vrc are the coupling matrix of the left L and 
right R leads with the central region. The main objective 
in the method is to calculate the Green’s function, which 
can be written as: 

G C {E) = (El - H c -Tj L — Si?)” 1 , 


(A2) 
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where / is the identity matrix, Tie = Vic 9i V^ c is the 
self-energy of each lead i = L,R, and gi = (E — H^) -1 
is the renormalized Green’s function of the semi-infinite 
lead i = L,R. In the linear response regime, the conduc¬ 
tance G is calculated as a function of the Fermi energy E , 
within the Landauer for mali sm. In terms of the Green’s 
function for the system Esmj q reac [ s: 


G(E) 


2e 2 

-T(£) 



T l G c T r G ] c , 


(A3) 


2 

V 



where T (E) is the transmission function across the con¬ 
ductor, and T £ = i[Ei — Ej] is the coupling between the 
conductor and the leads. 

Finally, within the Green’s functions formalism, the 
LDOS per atom is proportional to a Green’s function 
matrix element, 71 in a particular the LDOS at the atomic 
position j, can be expressed as Vj ( E ) = -4mG^ J ( E ), 
in such a way that the system LDOS is calculated by 
V (E) = -iTr pm Gc (£?)] = £ all atoms C7 


FIG. 9: (Color online) Normalized average LDOS per column 
of atoms as a function of the ^-position for the first two quan¬ 
tized /c^-values; (a) n = 1 at E = 0.048 70 and (b) n = 2 at 
E — 0.14 70 . Red and blue colors indicate the sublattice, and 
the orange lines indicate the position of the hydrogen atoms 
on the ribbon. 


Appendix B: Spatial Distribution of Localized and 
Resonant States 



FIG. 8 : (Color online) LDOS per atom, red and blue colors 
represent the sublattice of the atom. The bottom panel shows 
the LDOS (black line) and the conductance (red line) for a 
D2-configuration on a metallic AGNR N — 11 separated by 
L = 15. The arrows label the energies of the LDOS peaks 
plotted in the panels (a-d). 


We present the spatial distribution of the LDOS per 
atom where the red / blue color indicates different sublat¬ 
tice, for the lower sharp states belonging to the metallic 
band of the ribbon, labeled by (a-d) in the bottom panel 
of Fig. [8] The bar color scale represents the density 
probability per sublattice. It is possible to observe some 
general behavior for these states. First, there is an odd 
sequence of nodes along the longitudinal ribbon direction, 
defined by the allowed fc™-values. Second, the boundary 
conditions imposed by the impurities determine the spa¬ 
tial distribution of the LDOS for each sublattice, which 
exhibits the same feature of standing waves in a pipe with 
an open ends. This kind of behavior is similar to the ob¬ 
served in edge states of zigzag ribbons. To emphasizes 
these LDOS features, in Fig. [9] shows the normalized 
average LDOS per column of atoms (^/-direction) as a 
function of the x-position for the first two quantized k%- 
values over the envelope metallic band (corresponding to 
panels (a) and (b) of Fig. |8|. 

Figure [To] shows the spatial distribution of the LDOS 
per atom corresponding to the energy states labeled by 
(a-d) at the bottom panel of the figure. The (a) panel 
shows the mid-gap state located at zero energy, present 
in all ordered configurations, rega rdless if the ribbon is 
metallic or semiconductor. 27 39 40 This impurity induced 
state is highly localized in the nearest-neighbors sites to 
the adsorbed hydrogen atoms. 

The (b) panel shows a non-conductive interface state, 
which is determined by the band-folding in the /in¬ 
direction. This state corresponds to the first allowed 
discrete wavenumber, which is below the continuum of 
the graphene, and consequently, has a higher lifetime. 
This state is distributed between the hydrogen ad-atoms 
inside the doped regions. 

In panels (c) and (d) show two resonant and ex¬ 
tended states, determined by the band-folding in the k y - 
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direction, which are hybridized with the continuum of 
the graphene, shown in Fig. [4] The main feature of 
these states is that they are spread in a narrow energy 
window, so any electron from the leads i nject ed at these 
energies easily passes across the system.ESEH 
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FIG. 10: (Color online) LDOS per atom of: (a) the zero- 
energy state and (b-d) the interface states at high energies, 
for the same ribbon parameters of Fig. [8] Red and blue colors 
indicate the two graphene sublattices. 


Appendix C: Comparison with 
Density-Functional-Theory Calculations 


The DFT calculations are performed within the super¬ 
lattice approach using the density functional based tight- 
binding method implemented by the DFTB± code 73 -1with 
the associated Slater-Koster parameters.^ DFTB+ pro¬ 
vides an efficient quantum simulation tool based in DFT, 
which allows us to calculate, in affordable times, the elec¬ 
tronic properties of a system with the same size as cal¬ 
culated within a tight-binding approach. Self-consistent 
charge calculations are converged up to a tolerance of 
10 -8 e for the unit cell similar to the one shown in Fig 
TTJ Ribbons are repeated periodically using super-cell 


approximation and are separated by an empty space of 
10 A in the perpendicular directions. This cell is large 
enough to converge results with a single T-point. All 
atoms are relaxed within the conjugate gradient method 
until forces have been converged with a tolerance of 10 -3 
eV/A. We check that the relaxed geometry of a single 
hydrogen atom ads orbed on graphene reproduces well- 
known results G3E? 77 In agreement with these previous 
works, we found that the hydrogen atoms bonded to car¬ 
bon atoms are about 1 A over the graphene plane and its 
distance to nearest neighbor carbons is about 0.08 A. 


In Fig. 11 we show the isosurfaces of the second 


harmonic-like state calculated by using the tight-binding 
and the DFTB-f approach. The latter includes the geo¬ 
metrical corrugation induced by the sp 2 to sp 3 bonds of 
carbon atoms attached to hydrogen. These states, calcu¬ 
lated by two methods, have similar spatial distribution. 
In panel (a) we have plotted the spatial distribution of 
the LDOS per atom at E = ±0.46 70 ~ ± 1.2 eV obtained 
by the tight-binding method. As a comparison in Fig. |TT] 
(c) and (d), we have plotted the spatial distribution of the 
corresponding wavefunction (at energies E = 1.15 eV and 
E = —1.18 eV) obtained by using DFTB± calculations. 
To facilitate the comparison, atoms have been colored in 
black and red according the corresponding sublattice. It 
is important to mention that the energies of the latter 
states do not preserve the electron-hole symmetry due to 
the electronic correlation effects; however, they still have 
a good agreement with the ones previously obtained in 
the tight-binding model. 


Additionally, in the panels (c) and (d), a change in 
the sign of the wavefunction is evidenced around of the 
nodes. The wavefunction in (c) is even, while in (d) 
it is odd. The sign of resonant DFTB± wavefunctions 
changes from bonding-like to antibonding-like character 
for negative and positive energies from Ef , respectively. 
This parity change means that optical elements between 
both resonances will be certainly far from zero, and op¬ 
tical transitions between pairs of states could be also de¬ 
tected in experiments. 


The wave-functions arising from parabolic bands re¬ 
mind standing-waves on a rope, where both sublattices 
exhibit nodes at both ends, as shown in Fig. [lljb). This 
spatial distribution is different from the wave-function 
originated in the linear envelope-band of Fig. [9j where 
each sublattice presents a node at one end and a maxi¬ 
mum in the other. 


Finally, in DFTB± calculations we look at the energy 
separation between resonant states. To the separation, 
we associate an effective energy scaling y e ff, taking in to 
account the resonant levels far from the Fermi energy. We 
obtain a value of y e ff = 2.57 eV, in good an agreement 
with the one used through this work. 
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FIG. 11 : (Color online) Comparison between the LDOS per atom (red and blue colors represent the sublattice of the atom) 
calculated within the tight-binding model and the wavefunction at T-point calculated within the DFTB+ for a parabolic-band 
state. The (a) panel shows LDOS by the tight-binding model at E — ±0.46 70 ~ 1,2 eV, and in panel (b) the averaged LDOS 
over the y-axis as function of the ^-position. Red and blue colors indicate each graphene sublattice. The equivalent DFTB± 
wave-functions correspond to E = 1.15 eV and E = —1.18 eV in panel (c) and (d), respectively. 





















